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Abstract 

We review the status of QCD at hadron colliders with emphasis on precision predictions and the 
latest theoretical developments for cross sections calculations to higher orders. We include an 
overview of our current information on parton distributions and discuss various Standard Model 
reactions such as V7 ± /Z-boson, Higgs boson or top quark production. 



1 Introduction 



Historically, hadron colliders have explored elementary particle physics at the energy frontier the 
motivation being the discovery of new particles through direct production. This has been the case 
for the SppS at CERN leading to the discovery of the weak vector bosons as well as for Tevatron at 
Fermilab currently operating at a center-of-mass energy = 1.96 TeV with the discovery of the 
top quark. Shortly the Large Hadron Collider LHC at CERN with ^/S =14 TeV will commence 
operation which will realize a major leap forward in collision energy. Being long awaited the 
machine will allow access to the mechanism of electro-weak symmetry breaking, to search for 
the Higgs boson and, hopefully, it will open new avenues to test many proposed extensions of the 
Standard Model. To that end, two general purpose experiments ATLAS and CMS [1-4] as well 
as two specialized one, LHCb for 5-physics [5] and Alice for heavy-ion physics [6], have been 
installed. 
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Figure 1: Predictions for hard-scattering cross sections in pp collision at Tevatron and in pp collision at 
LHC as a function of the center-of-mass energy (WJ. Stirling in Ref. [7]). 



The expected cross sections for proton-proton scattering at LHC is large (see Fig.Q]). In particu- 
lar we will have large rates for many Standard Model processes such as the production of b-quarks, 
W ± and Z-bosons, jets (even with high cuts on the transverse momentum) and top quarks. Much 
of the physics is actually dominated by the gauge theory of the strong interactions, Quantum Chro- 
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modynamics (QCD). Comparing the rates for various processes at Tevatron and LHC in Fig.[H it 
is obvious that any search for physics beyond the Standard Model (BSM) like for superpartners in 
supersymmetric extensions (squarks, gluinos, ...), for Kaluza-Klein modes in models with extra di- 
mensions or even the search for the Higgs boson needs a very precise understanding of the known 
background from the Standard Model. Thus, new physics searches require precision predictions, 
most importantly in QCD. Moreover, the era of LHC implies a change of paradigm. We no longer 
test QCD, rather we use perturbative QCD as an essential and established part of our theory toolkit. 

Hard QCD is a large subject and, necessarily, the coverage here has to selective (see [7, 8] for 
other recent reviews on this topic). In this article, we will briefly review the physics concepts and 
theoretical framework for hard scattering reactions at hadron colliders, focusing on QCD. We point 
out achievements of the past years as well as open problems. We will briefly explain the property 
of factorization and discuss the parton luminosity in proton collisions. We will summarize the 
present knowledge on hard parton scattering cross sections for the production of W ± and Z gauge 
bosons, jets in QCD, heavy quarks, like top and bottom and the Higgs. We put particular emphasis 
on exact calculations of radiative corrections in QCD to next-to-leading order (NLO), next-to-next- 
to-leading order (NNLO) or beyond. Many other aspects such as e.g. the Monte Carlo approach to 
modeling hadronic interactions and parton showers we can only touch briefly or else, have to refer 
to the literature. 

With many Standard Model processes to be measured in the early days of LHC and the asso- 
ciated uncertainties to be understood on the way to the discovery of new physics [9], we hope that 
this review serves to illustrate a few aspects common to the underlying QCD dynamics. 



2 Perturbative QCD at colliders 

The basic prerequisite for the application of perturbative QCD at colliders is factorization for hard 
scattering processes. For hard hadron-hadron scattering this property implies that the constituent 
partons from each incoming hadron interact at short distance (i.e. at large momentum transfer Q 2 ). 
The property of QCD factorization rests on the fact that we can separate the sensitivity to dynamics 
from different scales (see e.g. [10]). Thus, for a cross section o pp ^x of some hadronic final state 
X in, say, proton-proton scattering we can write 

o pp ^x = / dxidx 2 dzfi(xui/)f j (x 2 ,iJ 2 )Oi j ^ k (xuX2,z,Q 2 ,a s (p 2 )^ 2 )D k ^x(z,i/), (1) 
ijk J 

where all functions have a clear physical interpretation. 

The parton distribution functions (PDFs) in the proton fi (z = q, q, g) describe the fraction x\ of 
the hadron momentum carried by the quark or gluon and the convolution of fi and fj determines 
the parton luminosity at the collider. The PDFs cannot be calculated in perturbation theory due to 
the proton being a very complicated multi-particle bound state. Rather, they have to be obtained 
from global fits to experimental data. The (hard) parton cross section &ij-+k depending on the par- 
ton types i, j and k is calculable perturbatively in QCD in powers of the strong coupling constant 
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Figure 2: Factorization for the hard-scattering cross sections in Eq. (Q]) in the QCD improved partem model. 



Oy and describes how the constituent partons from incoming protons interact at short distances of 
order o{\/Q). The final state X may denote hadrons, mesons, jets, etc. and needs another transi- 
tion from the perturbative hard partons in the final state to the observed particles. The necessary 
function D^x can therefore be a fragmentation function or also a jet algorithm. Here the interface 
with showering algorithms (based on a Monte Carlo approach) becomes particularly crucial. All 
quantities in Eq. © depend on the renormalization and factorization scale, fj r and /jf, which are 
usually taken to be the same. Throughout this review we set /j r = fjf = /j. The details of the inte- 
gration range in the convolution in Eq. (0Q) are controlled by the kinematics of the hard scattering 
process. Schematically QCD factorization can be depicted as in Fig.[2l 

Physical observables like the cross section o pp ^x in Eq. (OQ) cannot depend the factorization 
scale. In the perturbative approach, this implies that any dependence on fj in o pp ^x has to vanish 
at least to the order in a s considered. This property can be cast in the following form, 



din /j 



2°PP- 



X 



o(aj 



/+!- 



(2) 



It defines the commonly adopted approach to quantify uncertainties in theoretical predictions based 
on the variation of the renormalization and factorization scale. 

Let us briefly turn to hard scattering cross sections. There exist various approaches to the 
calculation of o pp ^,x ranging from easy to difficult as far as computational complexity is concerned 
as well as from inclusive to fully differential in terms of kinematical variables. First of all, there 
exist parton shower Monte Carlos (e.g. Herwig [11], Pythia [12, 13], Sherpa [14]) which are very 
important tools for understanding multi -parton scattering and the underlying event. 

For predictions building on exact matrix elements at leading order (LO), we have at our disposal 
many automated tree level calculations in the Standard Model, in its minimal supersymmetric ex- 
tension (MSSM) or in other BSM models utilizing programs like e.g. Alpgen [15], CompHEP [16], 
Helac-Phegas [17], MadGraph [18] or Whizard [19]. These tools provide first estimates for hard 
scattering cross sections through numerical phase space integration of the exact matrix elements. 
In this way, they are flexible as far as kinematics and the topology of a given hard scattering observ- 
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able is concerned and allow easy interfacing of LO calculations with parton shower Monte Carlos 
and, possibly, detector simulation. However, scattering reactions with exact matrix elements for 
more than 8 jets (particles) in the final state are currently at the edge of computational capabilites. 

At NLO level we do have some analytical (or numerical) calculations of Feynman diagrams 
yielding parton level Monte Carlos (e.g. NLOJET++ [20,21] or MCFM [22]). However, we have 
also seen recently significant progress based on string inspired techniques. At the edge of technical 
developments is the concept of exact NLO calculations interfaced with parton shower in programs 
as realized in MC@NLO [23,24], POWHEG [25,26] or VINCIA [27]. 

At higher orders in QCD perturbation theory, like NNLO some selected results are known 
mostly for inclusive kinematics. However, in view of LHC, we have witnessed significant progress 
in the last years to provide also predictions in completely differential kinematics [28,29]. Beyond 
this level of accuracy at, say, next-to-next-to-next-to-leading order (N 3 LO) only very few results 
are known, e.g. for deep-inelastic scattering (DIS) [30]. 

3 Parton luminosity at hadron colliders 

The parton luminosity in Eq. (QQ) is an indispensable ingredient of hard- scattering processes in- 
volving initial- state hadrons. At hadron colliders one has wide-band beams of quarks and gluons 
and, as is well known, the necessary PDFs of the proton fa (z = q, g) are not directly accessible 
in QCD perturbation theory. However, the scale dependence (evolution) of PDFs is governed by 
the splitting functions and predicted in a perturbative expansion in powers of a s . The universality 
allows for the determination of sets of PDFs in global fits to experimental data. Upon evolution 
this information from fits to reference processes can be used to provide cross section predictions 
at LHC energies and we can quantify the present uncertainties. 

3.1 Parton evolution 

The parton distributions in the hadron are distinguished by the flavor quantum numbers, which are 
additive. The valence distribution originates from differences of quarks and anti-quarks q — q. The 
proton is composed of the sea distribution (i.e. the sum over all flavors q + q) and of the gluon g. 

The independence of any physical observable on the scale /j immediately gives rise to evolution 
equations for the PDFs ft, i = q,q,g. From Eq. © we find that the scale dependence of /, is 
governed by 

d ( U(x,v 2 ) \_ V fdzf p mM) ^ ( W^M 2 ) \ ™ 

dlnS V fg&f) J yJ z\ P gqi (z) P g8 (z) J \ f 8 {x/z^ 2 ) J ' { ) 

J X 

which is a system of coupled integro-differential equations corresponding to the different possible 
parton splittings, see e.g. Fig. [3] where some Feynman diagrams contributing in leading order 
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Figure 3: Sample of Feynman diagrams for parton-parton splitting in leading order QCD. We indicate the 
collinear momentum flow (p incoming and xp outgoing) as it enters the calculation of the corresponding 
splitting function P^. See e.g. Ref. [31] for an operator definition of parton distributions. 
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Table 1: The number of Feynman diagrams contributing to parton (q,g)-boson DIS (vector bosons y/W 
or scalar (()) up to three loops. The NNLO splitting functions Py have been determined from the collinear 
singularity of these scattering reactions (see [32, 33]). 

QCD are displayed. The splitting functions Py, i.e. the kernels of these differential equations are 
universal quantities and can be calculated in perturbation theory from the collinear singularity of 
any hard scattering process. Thus, P has an expansion in powers of a? as 

p = a s P {0) +a 2 P w + a 3 s P {2) + ... , (4) 

where we have suppressed parton indices. The first two terms are needed for NLO predictions, 
which is the standard approximation, although often still with large uncertainties. Currently, the 
splitting functions are known to NNLO and in Tab. Q] we give the number of Feynman diagrams for 
the corresponding hard parton reactions in DIS from which the NNLO expressions P^ in Eq. © 
have been calculated [32, 33]. 

Physically, the evolution Eq. © states that one becomes sensitive to lower momentum partons 
as the resolution of the proton is increased, i.e. as the scale /j becomes larger. Given an input 
distribution at a low scale, say Q 2 = 10 Gev 2 , which has to be determined in a global fit from 
comparison to data, one can solve Eq. © to predict the PDFs at a high scale (see Fig.|4j>. Solutions 
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of Eq. © can be obtained by a variety of methods with available codes [34-36] and benchmarks 
are provided in Refs. [37, 38]. 
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Figure 4: Evolution of the valence, sea and gluon momentum distributions xf(x,Q) in the proton from 
a low scale at Q 2 = 10 GeV (left) to LHC energies at Q 2 = 10 4 GeV (right) for the parameterization of 
Ref. [39] showing the strong rise of the gluon at small x. 



Modern parameterizations of parton distribution from global fits account in particular for the 
effects of experimental errors and come with the according uncertainties, see e.g. the frame- 
work LHAPDF accord [37,40,41]. Much of the needed experimental information originates from 
deep-inelastic scattering data on structure functions from HERA for e ± p- scattering (HI, ZEUS) 
and from fixed targets (proton and deuterium) for up and [id scattering (BCDMS, NMC, SLAC, 
E665), as well as (anti-)neutrino-proton scattering (CCFR), see e.g. [42]. These data determine the 
quark distributions for light flavors at all momentum fractions x and through the scale evolution in 
perturbative QCD the gluon distribution at medium and small x. Further information on the flavor 
content of the nucleon is provided by structure function data for F£ m from HERA for the charm 
distribution and by Drell-Yan data on proton-nucleon targets (E605, E772, E866), which deter- 
mine the sea quark distributions, in particular fa and fj. The Tevatron experiments CDF and DO 
are able to constrain the ratio f u / fd at high x with the rapidity asymmetry in W-boson production 
and the gluon distribution at high x with the help of inclusive jet data. More recently, by relaxing 
the assumption f s = fs, also information on the strange asymmetry f s and has been extracted 
from v(v)p scattering (NuTeV, CCFR). 

3.2 Parton distributions from HERA to LHC 

Given the great importance of deep-inelastic scattering data a question that has been frequently 
addressed in the past is, of course, the impact HERA data for LHC predictions, in particular as far 
as the parton luminosity is concerned (see Refs. [38,43] for further discussion). 

As illustrated by the compilation of measurements in Fig. [5] (left), the data on the structure 
function F2 in deep-inelastic e ± p scattering extend over a wide range in x and Q 2 . Considering, on 
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Figure 5: Left: The structure function F2 in deep-inelastic e p scattering as a function of x and Q 2 (from 
Ref. [44]). Right: Parton kinematics at LHC and at HERA (W.J. Stirling in Ref. [7]). 



the other hand, the allowed region for parton kinematics at LHC in Fig. [5] (right), it is clear, that 
there is a large overlap in x with the range covered by HERA. However the relevant hard scale Q is 
typically two to three orders higher due to the increased center-of-mass energy \/S. Details of the 
parton kinematics at LHC depend, of course, on the invariant mass M of the final state and on the 
rapidity y. The dominant values of the momentum fractions are xi.2 ~ (Me ±y ) / y/S and variation 
of M and y at fixed y/S tests the sensitivity to partons with different momentum fractions. 

The large difference in the hard momentum scale Q between HERA and LHC requires the 
parton evolution based on Eq. © to be sufficiently accurate in perturbative QCD. The necessary 
perturbative accuracy for quantitative predictions is approached at NNLO [32,33]. The stability of 
evolution is shown in Fig.[6l where the scale derivatives of quark and gluon distributions at fi 2 ps 30 
GeV 2 are displayed. Obviously, the expansion is very stable except for very small momentum 
fractions x < 10~ 4 which shows that the perturbative evolution Eq. ([3]) is applicable down to very 
small x. In terms of LHC parton kinematics, this corresponds to perturbative stability for central 
rapidities \y\ < 2, while modifications are at most expected in the very forward (backward) regions 

|y|£4. 
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3.3 W and Z-boson production at LHC 

The immediate question arises: What is the impact of our current knowledge of parton distribu- 
tions on the precision of LHC predictions, for instance for W ± , Z-boson rapidity distributions, 
which often have been considered "standard candle" processes for the parton luminosity [38,45]. 
The corresponding cross sections are known to NNLO in perturbative QCD, and according to 
Eq. © one can quantify the theoretical uncertainties obtained by varying the renormalization and 
factorization scale ^ by the conventional (although arbitrary) factor of two around My/ z- The per- 
turbative stability of the results in Fig. [7] nicely demonstrates the necessity of considering higher 
order perturbative corrections through NNLO in QCD. It would be impossible to make precision 
predictions, or perform precision analyses, based solely on the rough (and non-overlapping) LO 
and NLO error estimates. 

Recent improvements in the parameterizations of PDFs, though, have been shown to signif- 
icantly affect predictions for physical cross sections at LHC. An independent treatment of the 
strange quark distributions f s and fg (hence their uncertainties), for instance, has an impact on the 
correlated uncertainties of the light sea quarks, because neutral current deep-inelastic data on F2 
constrains the combination 4/9 (/» + /») + 1 /9 (fd + fj + f s + fs) ■ In consequence, the size of the 
uncertainty on the sea quarks for values x ~ 10~ 3 — 10~ 2 at hard scales Q 2 ~ roughly doubles 
from ~ 1.5% to ~ 3% for the MSTW group [50, 51]. CTEQ in their newer sets (e.g. CTEQ6.6) 
has improved the treatment of the charm contribution to the deep-inelastic structure function F2 at 
HERA by implementing now a general-mass formalism for a variable flavor number scheme con- 
sistent with QCD factorization, see [52, 53]. The reduced charm component of F2 is compensated 
by larger light quark distributions f u and fy at small x as illustrated in Fig.[8l 
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Figure 7: The rapidity-dependent cross sections for gauge-boson production at the LHC, using the partons 
of Ref. [46] and estimates of the theoretical uncertainty from variations of the scale /u (from Ref. [47]). 



As an upshot, the predictions for W - and Z-production cross sections at LHC being sensitive 
to PDFs in the x ~ 10 3 range shift by 8% between the sets CTEQ6.6 [48] and CTEQ6.1M [49]. 
Although this particular shift originates from theoretical improvements long overdue, it is an ex- 
ample that PDFs and their associated uncertainties will have a significant impact on the precision 
of the 'gold-plated' W ± - and Z-cross-sections and W ± /Z-ratio calibration measurement. In this 
context, it should also be stressed, that PDF uncertainties in the region of very small momentum 
fractions, x ~ 10~ 5 , (as e.g. displayed in Fig. [8]) largely rely upon extrapolations of data and 
represent a certain parameterization bias. 



3.4 Parton distributions and the search for new physics 

Apart from gauge boson production, there are prominent measurements at LHC which depend 
on our knowledge of parton distributions and, in turn, might be used to improve it. High-EV jet 
cross-sections, for instance, are a particularly prominent place to look for BSM effects. The dis- 
covery of new physics, such as e.g. jet signals for low mass strings [54], large extra dimensions 
or models parameterized in terms of contact interactions becomes sensitive to the uncertainty of 
the gluon PDF especially at low-x. Recently, also top-pair-production has been proposed as an ad- 
ditional calibration process at LHC, because its PDF dependence is anti-correlated with Z-boson 
production [48] and correlated with Higgs boson production, especially for larger Higgs masses. 
Presently, however, the sizable theoretical uncertainties at NLO in QCD are limiting the applica- 
bility of this proposal. 

For di-jet rates at LHC, e.g. the consequences for predictions from large extra dimensions have 
been analyzed [55] through a modified renormalization group equation for the strong coupling oc 5 , 
where the running of 0^ accelerates due to power corrections as the compactification scale M c of 
the large extra dimensions is approached. Results of a study for the di-jet transverse momentum 
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Figure 8: The distributions /„, f c / and f g of the CTEQ6.6 fit [48]) with the estimated error bands (shaded 
area) at the scale /u = 100 GeV. Comparison with previous results (CTEQ6.1M) [49] is denoted by solid 
lines (from Ref. [48]). 



distribution (p t ) with the event generator Pythia are displayed in Fig. [9j The plot clearly indi- 
cates the reduced sensitivity to extra-dimensions because of the PDF uncertainties in the di-jet 
/^-distribution. Hence there is need to either look at ratios of rates, <s pp ^^ j e t s over o pp ^2 jets or ' 
else at angular correlations of di-jets to reduce the parton luminosity dependence. 

On the other hand, PDF uncertainties most likely do not affect the discovery potential of a 
Higgs in the mass range 100 — 1000 GeV or a high mass Z' in the mass range 150 — 2500 GeV. 
Apart form the hadronic di-jets, promising other measurements to be conducted at LHC itself 
also include direct photon production to constrain the gluon PDF at low-x or the V7 ± -asymmetry 
to obtain information on the low-x valence PDFs. In particular, the improved description of the 
(anti-)strange quark distributions leads to interesting implications for collider phenomenology. For 
instance the production of a charged Higgs boson H + via the partonic process c + s^H + , provides 
an example of a BSM process that is sensitive to the strange PDF in models with two or more Higgs 
doublets. The cross section also depends on a possible intrinsic charm component of the proton 
and the recent PDF set CTEQ6.5c provides various models for such a component [56]. 

4 Parton cross sections 

Let us next turn to the hard scattering cross sections initiated through the constituent partons j 
of the incoming protons. In general, we consider &ij->x where X denotes any final state allowed 
by the Standard Model or its possible extensions. Calculations of <Jj/-+x result in predictions for 
experimental signatures and, eventually, determine the power to discriminate BSM signals or, e.g. 
Higgs boson production, from known background of the Standard Model. The key issue for QCD 
theory is the reliability of signal and background estimates. 
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Figure 9: Di-jets cross sections at LHC as a function of p t for two compactification scales M c = 2 TeV 
(left) and M c = 4 TeV (right) of the extra-dimensions. Shown are different predictions for various numbers 
of extra dimensions, i.e. 2,4 and 6, and the Standard Model zone incorporates the PDF uncertainties of 
the CTEQ6.1M set. The horizontal line shows the sensitivity limit corresponding to an LHC luminosity of 
100 fbr 1 (from Ref. [55]). 



4.1 QCD @ NLO 

At a hadron collider, the problem of signal significance has various aspects. Experimentally, mea- 
surements of hard scattering reactions require reliable identification of leptons (electrons, muons), 
a good understanding hard jets at high transverse momentum (p t ), especially £-quark jets, and a 
sufficiently precise calibration of the jet energy scale. Moreover, BSM or Higgs searches rely heav- 
ily on the presence of large missing transverse energy {f t ) to reject Standard Model background 
compared to the signal. 

On the theory side, as briefly mentioned in Sec. [2l we have various levels of accuracy for the 
hard scattering process (assuming that the underlying event and multiple parton interaction are 
modeled by shower Monte Carlos). Estimates to LO in QCD based on exact matrix elements seem 
mandatory in search scenarios for studies of distributions, e.g. in p t or the (pseudo-)rapidity (r|) 
and for assessing the effects of kinematical cuts. It is well known, that the overall normalization 
and, in particular, the hard tail of these distributions (e.g at high p t or f t ) are not well modeled 
by shower Monte Carlos alone. However, any LO prediction has large theoretical uncertainties, 
typically estimated by the scale variation, Eq. ©. Consider, for instance, the cross section for 
pp — > W jets, which is of O(aj) at LO. From a variation of the coupling of A(oc^°) ~ 10% one 
can roughly estimate a cross section uncertainty of A(o LO ) ~ 40%. Thus, one needs to go beyond 
the Born approximation for certain processes. 
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Table 2: Scattering processes at LHC for which the radiative corrections to NLO in QCD are needed, as 
summarized in Les Houches 2005 (from Ref. [62]). 



Perturbative QCD corrections at NLO to scattering processes are essential for the rates and 
shapes of distributions of Standard Model processes as well as for BSM searches, where they may 
have an impact on the signal significance. Often, one encounters large K-factors and also new 
parton channels open up at NLO which may eventually dominate beyond tree level, a prominent 
example being single-top production (see Sec. 15.21) . In a series of workshops a number of key 
processes at LHC has been identified which need to be known to NLO in QCD. These are sum- 
marized in the so-called LHC "priority" wishlist in Tab. |2] and the computation of these radiative 
corrections is presently a very active field of research. 

real corrections — subtractions virtual corrections + subtractions 

(IR divergent) (IR divergent) 

cancellation of singularities 

finite partonic cross sections 

phase space integrations 

convolutions with PDFs 
Monte Carlo 

Figure 10: Outline of a generic calculation of NLO QCD corrections to a (mufti -particle) scattering process 
in the traditional approach. 



The obvious question is, of course, why the calculation of one-loop corrections in QCD is so 
difficuft? After alf, the conceptuaf issues are aif sofved and any computation can folfow a straight- 
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forward algorithm: Draw all Feynman diagrams and evaluate them, then use standard reduction 
techniques for tree and loop amplitudes. While this is true in principle, it is hard in practice with 
known bottlenecks, because intermediate expressions are much more complicated than the final 
result. Thus, let us look at the outline of a generic NLO calculation as displayed in Fig. [TOl For 
the scattering reaction 2 — > n partons, the basic ingredients in the cross section calculation are 
the real corrections, i.e. the tree level 2 — > (n+1) parton reaction do real , and the one-loop vir- 
tual corrections to the 2 — > n parton amplitude, which are subject to the standard ultraviolet (UV) 
renormalization. 

The latter contribution to the cross section, that is the one-loop virtual correction Jo virtual , 
in the standard Feynman diagram approach generates large expressions, although we do expect 
large cancellations between the diagrams in a gauge theory as a consequence of gauge invariance. 
Specifically, one is required to calculate tensor integrals like, e.g. 



d D 



Pi 



(pj-mlXfa-ktf-ml) 



(5) 



Unfortunately, these become rather complicated for five or more external particles, one well-known 
problem being the numerical stability for all allowed configurations of the external momenta. 
Without going into details here, suffice it to say, that reduction algorithms for tensor integrals 
are some 30 years after the work of Passarino and Veltman [63] still a very active field of re- 
search, see e.g. [64-66]. Moreover, no completely general libraries are available here (see e.g. 
LoopsTools [67] for public code). 

Coming back to Fig. \W[ a characteristic feature of both contributions <ia real and Ja virtual is 
the presence of infrared (IR) divergencies due to soft and collinear regions in phase space. The 
physical cross section o^9 being the sum of both is, of course, IR finite after absorbing the 
initial state collinear singularities into the PDFs by mass factorization (see e.g. [68] for details). 
As we aim at parton level Monte Carlos to NLO accuracy with a flexible phase space integration 
allowing for kinematical cuts, the IR divergencies need to be treated accordingly. Among the many 
proposed methods (see e.g. [8]), the so called dipole subtraction [69,70] has emerged as a standard 
procedure. Here, the cancellation of infrared singularities due to collinear partons or soft gluons 
is implemented "locally" by subtracting (over the entire phase space) functions (the dipoles) that 
approximate the singularities of the real emission part. Subsequently, the integrated dipoles are 
added to the virtual corrections. Employing dimensional regularization ( d = 4 — 2e), the master 
formula reads [69,70], 



_NLO 



n+1 



do 



real 



e=0 



Jo dipole 



)e=0. 




^virtual + J Ja dipole 









(6) 



Je=0 



which is understood to contain also the mass factorization of the remaining initial state collinear 
singularities. With Eq. © one arrives at a finite partonic cross section o^9 which can be inte- 
grated numerically over the available phase space and convoluted with the PDFs. Let us stress 
however, that for any practical solution of Eq. © speed and stability of the numerics are criteria of 
paramount importance. Quite often for instance, the programs for a particular scattering process 
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need dedicated optimization. Thus, presently there is a lot of room for technological progress with 
respect to automatization and algorithms in order to fill the empty spaces in Tab. |2l 

4.2 New theory developments 

One specific direction of research in theory during the past few years has been towards new ana- 
lytic techniques to calculate gauge theory amplitudes. The specific focus has been on a recursive 
approach in which all intermediate quantities are on-shell and hence gauge invariant (see [71] for 
a recent review). As a matter of fact, techniques for computing tree amplitudes recursively are 
well established since a number of years [72]. Moreover, it has been realized that an efficient 
management of the quantum numbers for a given scattering amplitude reduces the computational 
complexity by far. The known methods include so-called color ordering, the use of helicity ampli- 
tudes and a decomposition of QCD amplitudes exploiting effective supersymmetry (SUSY) (see 
e.g. [73]). In addition, factorization properties of amplitudes in the soft and collinear limits serve 
as a strong check (see for instance [69,73]). 

In a helicity basis, amplitudes for scattering processes are classified according to the number 
of '±'-states of the external partons. The so-called maximal helicity violating (MHV) amplitude 
denotes the configuration with the largest difference of '+' and '—'-states, e.g. n — 2 for a n- 
gluon amplitude A„ at tree level. The tree level (color ordered) n-gluon MHV amplitude Aj ? ree 
takes a particularly simple and elegant form [74]. In terms of (Weyl) spinor inner-products (jl) = 
uZ{kj)u + {ki) for massless Weyl spinors u±(k) of momentum k we can write (see also e.g. [75]), 

A ^ 1 " 2 "' 3+ » + ) =i (12)ff.l.W (?) 
which is in a certain sense "all-order" information, because Eq. © holds for any number n of 
external gluons. 

The importance of helicity amplitudes became again apparent upon applying twistor space 
methods [76] and by realizing that tree amplitudes A* ? ree in gauge theories possess unique analytic 
properties which become manifest, if considered as functions of complex momenta k, i.e. under 
a shift k — > k(z) for a complex valued parameter z- These analyticity properties of Aj[ ee (z) can be 
turned into recursion relations for the case of gluons [77,78] as well as quarks and scalars [71, 
79, 80], which allow the construction of rc-point helicity amplitudes from on-shell (n — 2)-point 
amplitudes. From Aj[ ee (z) the physical amplitude at z = 0, i.e. A{ 3 ree (0) can be reconstructed simply 
by exploiting Cauchy's theorem of complex analysis. The key feature of complex kinematics is 
the fact that the three-parton primitive MHV amplitude, e.g. Eq. © for n = 3, does not vanish 
for on-shell complex momenta k(z), while it does so for real k. For details of the complex shift 
k — > k(z) we refer to [78]. 

The real quest, of course, has been in devising improved methods for the calculation of one- 
loop corrections to scattering amplitudes. This has been a very active field of research over the 
past few years and it has been realized that unitarity methods provide the additional key ingredi- 
ent here. Unitarity is a fundamental concept in quantum field theory which manifests itself for 
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Figure 1 1 : Schematic representation of the unitarity approach: Sewing of tree level amplitudes leads to the 
imaginary part of one-loop n-point amplitudes. The real part is determined subsequently from the analytic 
properties of the amplitudes (adapted from Ref. [81]). 



instance in the Cutkosky cutting rules for Feynman diagrams (related to the optical theorem). In 
the computation of one-loop amplitudes, unitarity appears as fusing rules for amplitudes [71, 82]. 
One aims at reconstructing the real part of a given one-loop amplitude from the imaginary one by 
sewing together tree level amplitudes. To that end, one uses the fact that any one-loop amplitude 
can be expressed in a basis of scalar integral functions, i.e. boxes, triangles, and bubbles. This 
is obvious from the standard reduction techniques, e.g. [63] and it is sketched schematically in 
Fig. QT] Unitarity cuts then allow to identify uniquely the contribution to the individual integral 
(i.e. the coefficients c^j, C3 j, C2 : k m Fig- ECD from its imaginary part, for example a triangle loop 
from calculating in\n(— s) — > In 2 (—s). 

Generalized unitarity (see e.g. [83]) developed further the idea of reconstructing the coeffi- 
cients C4 u etc. by imposing quadruple cuts, which constrain all components of a given one-loop 
integral ©. As an upshot, all terms with logarithmic dependence are cut-constructible. One-loop 
QCD amplitudes, however, also contain rational (non-logarithmic) terms, which are rather difficult 
to derive and require substantially more effort. 

To that end, as a further principle to organize the calculation a SUSY inspired decomposition 
of one-loop amplitudes has been very useful. For the amplitude A„ with n external gluons and a 
gluon circulating in the loop we can write down, 

Al = (Ag + 4^[ + 3Aj;) - 4 (a{+A*) + A s n , (8) 

5\£ = 4 SUSY 91 = 1 chiral SUSY 9i = scalar 

where the internal loop degrees of freedom are ordered in terms supersymmetric multiplets, i.e. 
the fA£ = 4 multiplet (1 gluon, 4 Weyl fermions, 6 real scalars), the chiral fA£ = 1 multiplet (1 
Weyl fermion, 2 real scalars) and the 2\£ = part with a complex scalar in the loop. This has the 
great advantage that 9£ = 4 and the 9l = 1 contributions are completely cut-constructible. For the 
rational part in the 3\£ = bit, more involved recursions have been developed. 

Finally, the computation of the one-loop virtual corrections to the six gluon amplitude has been 
completed through the effort of many groups [84-95] (see Tab. [3]) and the analytic results have been 
confirmed for specific phase space points by a completely numerical evaluation [96] . Recently, the 
numerical approach has been developed further and shown to have promising potential [97]. We 
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cut 
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- + + ++ 


[84] 
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[90] 
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[84] 


[85] 
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[84] 


[85] 


[87] 


[94,95] 


- + ++ 


[85] 


[86] 
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[93] 


- + -++ 


[85] 


[88,89] 


[92] 


[95] 


-+-+-+ 


[85] 


[88,89] 


[92] 


[95] 



Table 3: The analytic analytic computation of the one-loop QCD corrections to the individual helicity 
configurations of the six-gluon amplitude in the decomposition of Eq. ([8]) as a community effort (adapted 
from Ref. [81]). 

therefore expect that the NLO correction to the four-jet cross section at LHC are within sight. 

Let us end the discussion by mentioning a few directions for further development. Clearly, the 
new techniques have to be employed in complete cross section calculations following the steps 
outlined in Sec. 14.11 Moreover, at LHC many processes of interest contain either gauge boson 
or bottom and top quarks (see Tab. |2]). Thus, the formalism sketched above needs to be carried 
over to case of massive (colored) particles. This requires a number of extensions, be it the helicity 
formalism or the methods for calculating massive one-loop integrals from generalized unitarity 
cuts (see e.g. [98,99]). We should also mention, that of course many other developments in theory 
have pushed the precision frontier for QCD predictions further, be it for multi-loop calculations 
(NNLO and beyond), for resummations or simply for algorithms and tools. Unfortunately we 
could not touch those aspects. 

4.3 Complete NLO and NNLO results 

Let us conclude this Section by summarizing the state-of-the art for complete NLO cross section 
predictions for many-particle production. It has become clear in the preceeding discussion (see 
Tab. [2]) that current theory research is focused on processes with 2^3 and 2^4 potentially 
massive particles like top quarks, V7/Z-bosons, etc.. 

A few outstanding results that have appeared in the last one or two years in this respect ad- 
dress for instance Higgs production in the Standard Model. Here a number of reactions has been 
investigated at NLO, such as pp — > H + 2 jets via gluon fusion [59] (in the heavy top limit, see 
Sec. |6]) and via weak interactions [100] as well as pp — > H + 3 jets in vector-boson fusion [101]. 
These results provide important information e.g. for the extraction of the Higgs coupling to vector 
bosons at LHC. The production of vector bosons has been considered to NLO accuracy in QCD, 
e.g. for the reaction pp — »■ VV + 2 jets via vector boson fusion [102-104], for tri-boson production 
(pp -> ZZZ.WWZ) [60, 61] and for pp -> WW + 1 jet [57, 58]. Especially, the latter case is an 
important background for Higgs production in the low mass range (H —> WW) and subsequent 
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semi-leptonic decay of the W -bosons. A largely complete list of NLO calculations for LHC pro- 
cesses including electroweak corrections and also a number of new physics signals with e.g. the 
SUSY QCD corrections has been given in Ref. [62]. 

Some predictions for pp — > 3 particles are published as public codes. Many of these programs 
also provide continuous updates which are made available to the community e.g. through the 
CEDAR project [105]. To mention a few explicitly, let us point out NLOJET++ [20,21] which 
comes as a multipurpose C++ library for calculating jet cross sections at NLO, e.g. for three-jet 
rates in hadron collisions. The program MCFM [22] calculates for instance the production of 
a gauge boson or a Higgs in association with jets at hadron colliders. The PHOX family with 
DIPHOX [106] and JETPHOX [107] deals specifically with hard QCD radiation of photons along 
with jets in hadron collisions, which is another important background for the low mass Higgs, e.g. 
in the di-photon mode (H — > yy). Predictions for the photon pair- production background including 
a resummation of the transverse momentum p t of the di-photon pair [108, 109] have been subject 
of the most recent improvements in the program ResBos [110]. 




Figure 12: Sample of Feynman diagram (hexagon with internal masses) for electroweak corrections to 
e + e~ — > 4 fermions. 

NLO radiative corrections to scattering reactions with six external particles, i.e. 2 — > 4 pro- 
cesses, constitute the current technological frontier. Benchmark results are the calculation of the 
complete electroweak corrections to — > 4 fermions [111,112]. This reaction involves at the 
loop level extremely difficult hexagon integrals with masses, see Fig. [[21 Of similar complex- 
ity are the NLO electroweak corrections to Higgs production in association with a neutrino-pair, 
e+e — >■ vvHH, obtained by the GRACE group [113, 1 14]. Also the NLO QCD corrections to the 
combined production of a top- and a bottom-pair, i.e. the process yy — > ttbb, are known [115]. 

Finally, there is of course demand for fully differential QCD predictions to NNLO for hadron 
collider processes. As mentioned above, this scope has been achieved e.g. for the di-lepton pair 
production in Drell-Yan [28] or Higgs production in gluon fusion [29] together with the parton 
evolution [32, 33]. However, it remains a challenge for hadronic di-jet production, where large 
statistics even with high-p f cuts is anticipated at LHC. The measurement of gluon jets would 
constrain for instance the gluon PDF at medium and large x and di-jet angular correlations are 
important observables in BSM searches for quark sub-structure. NNLO predictions for di-jets 
are likely to reduce the scale uncertainty and to improve the modeling of jets. Recent extensions 
of Tab. [2] also list the NNLO corrections to Higgs production in vector boson fusion, to top-pair 
production and to V + ljet, where V e {y, W ± ,Z}. 
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Unfortunately, the calculation of NNLO cross sections is very difficult. Although many (two- 
loop) virtual amplitudes are known since some years, the cancellation of IR divergencies between 
virtual and real corrections remains highly non-trivial and the numerical phase space integration 
very difficult, i.e. the NNLO equivalent of steps outlined in Fig. \10\ This is a vast subject on its 
own, which we will not pursue further here (see e.g. [8] for a brief review). Suffice it to say, that 
progress in this direction has been achieved only recently for the differential distributions in the 
case of e + e~ — > 3 jets [116, 117]. 

5 Top quark production at LHC 

Top quarks will be copiously produced at LHC. For the pair-production mode the collider will 
accumulate very high statistics of approximately 8 • 10 6 events with /f-pairs with 10 fb _1 per year 
in the initial low luminosity run [2, 3]. This data will allow for numerous measurements, e.g. of 
the top-mass, where the experiments aim at an accuracy of Am t = o(l)GeV, and also for tests of 
the production and the subsequent decay mechanism including anomalous couplings and top-spin 
correlations (see e.g. Ref. [118]). Top quark decay (t — > Wb) leads to very characteristic signatures 
allowing for event reconstruction in many channels through the observed leptons, (^-flavored) jets 
and missing f t . In general, top quarks make up large part of background for Higgs production or 
BSM searches. Moreover, due to the large mass (currently m t = 170.9 ± 1 . 1 (stat) + 1.5 (syst) GeV, 
see [119]) close to the scale of electroweak symmetry breaking top quarks play a prominent role 
in many new physics scenarios (see e.g. [120]). 

5.1 Top quark pair-production 

The hadronic heavy quark pair-production is known to NLO in QCD for many years [121, 122] and 
it still serves as an example to illustrate many generic features of QCD corrections to the production 
of heavy colored particles, see Fig. [13] for the corresponding Feynman diagrams to leading order. 
Depending on the collider, i.e. Tevatron (pp) or LHC (pp), the parton luminosities enhance the 
respective parton channels qq and gg. Thus, at Tevatron, ^-annihilation saturates the total cross 
section to (90%) , while at LHC a similar dominance of gluon-fusion holds. The channels qg(qg) 
newly opening up at NLO contribute only 0(1%) at the scale /j = m t at both colliders. 




Figure 13: Complete set of Feynman diagrams to leading order for the heavy quark pair-production in 
light quark annihilation, q + q — > t + 1, (left) and in gluon fusion g + g — > t + 1 (three diagrams on the right). 
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Figure 14: The tt total cross section to NLO QCD as a function of m t for Tevatron at y/S = 1.96 TeV and 
CDF data [123] for m t = 171 GeV (left) and LHC at y/S = 14 TeV (right). The solid line is the central value 
for /u = m t , the dashed lower and upper lines correspond to fj = 2m t and n = m t /2, respectively. The band 
denotes the additional PDF uncertainty of the CTEQ6.5 set [124]. 



In Fig. [[4] we show the NLO QCD predictions for the total cross-section of ff-pair production 
at Tevatron and LHC. The band denotes the scale variation in the usual range (m t /2 > /j > 2m t ) 
as well as uncertainties related to the parton luminosity. At Tevatron, the theory error budget at 
NLO in QCD is slightly asymmetric +12%/ — 15% which breaks down to a scale uncertainty of 
+5%/ - 10% and a PDF uncertainty of +7%/ - 5%. The latter one is due to the PDFs being 
sampled in the large-* region, where especially the gluon is poorly constrained. At LHC, the 
total theory error is 15% which consists of a scale uncertainty of 11% and a much smaller PDF 
uncertainty of 4%. Here, the cross section is sensitive to the gluon PDF in a range well covered 
by HERA (see Fig. [5]). Different sets of global PDFs agree within the given error bands, although 
it should be pointed out that there can be sizable shifts in the central values. For example, there is 
a 3% shift in the central value between the CTEQ6.5 and CTEQ6.6 sets [48] with correct heavy 
flavor treatment and the older set CTEQ6.1M (see also [125] for a recent discussion). 

The theoretical prediction can be improved in specific kinematical regions. Near threshold a 
Sudakov resummation can be performed [126-128], which stabilizes perturbative predictions if 
the /f-pairs are produced close to partonic threshold as for instance at Tevatron and, perhaps, to 
a lesser extent at LHC. Further improvements of the theoretical accuracy need the NNLO QCD 
corrections, which are mandatory for a precision of better than 0(10%) as envisaged by the LHC 
experiments. First steps in this direction have been undertaken by evaluating the interference of the 
one-loop QCD corrections [129] and by deriving the virtual contributions to heavy-quark hadro- 
production at two loops in the ultra-relativistic limit m 2 <C [130, 131] based on a simple 
relation of massive and massless amplitudes in the limit m — > [132] (see also the review [133]). 
A precise understanding of the kinematical region m — > beyond NLO is of immediate relevance 
also bottom-pair production over a large kinematical range and heavy flavor production at large p t . 
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Top quark mass determinations at LHC are usually planned to proceed through direct recon- 
struction because the theoretical accuracy of the total cross section is presently insufficient, see 
Fig. [H An interesting alternative (see e.g. [3,4]) involves a mass measurement through J/\\f final 
states from ^-decays in the top quark decay chain (t — > Wb). The y/y-reconstruction is supposed to 
give an accurate measurement of the Z?-quark momentum thanks to the relatively high mass of the 
meson. However, currently, it is also limited by our knowledge of the heavy-quark fragmentation 
the perturbative description of which has been extended to NNLO in QCD only recently [134-137]. 
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Figure 15: Left: The tt invariant mass spectrum at LHC for m t =170 GeV at NLO in QCD together with 
the scale (dashed) and the PDF (dotted) uncertainties for the CTEQ6.1M set. Also plotted are predictions 
at LO in QCD (normalized to the NLO total cross section) with (dark dash-dotted) and without (light dash- 
dotted) the NLO electroweak corrections for the CTEQ6L1 set (from Ref. [120]). Right: The tt invariant 
mass spectrum at LHC including ^-channel graviton exchange and the effect of a couple of Kaluza-Klein 
resonances in an extra dimensions model (from Ref. [120], where the model parameters are specified). 

Other observables of experimental interest at LHC are differential distributions in the transverse 
momentum p t of an identified top quark, or in the invariant mass m t j of the ff-pair. (see e.g. [120]). 
Predictions for the latter are displayed in Fig. [15] (left) based on NLO in QCD as calculated e.g. 
with MCFM [22] and also including NLO electroweak corrections [138] for the CTEQ6.1 set [49]. 
In the TeV-region for the invariant mass m t j both the NLO QCD corrections and also the elec- 
troweak radiative effects grow. The dominant theoretical errors however come from the scale and 
the PDF uncertainties. In particular the latter start to increase because the dominant contributions 
come again from the poorly known large-* region. 

It has been pointed out [120] that the m f f-distribution also provides a window to new physics, 
where s-channel resonances may become visible. Fig. [15] (right) nicely illustrates the effect of 
graviton exchange in a model with one extra dimension compactified to a S 1 /Z2 orbifold (see 
e.g. [139] and references therein). Such a model leads to a tower of Kaluza-Klein modes giving 
rise to a series of resonances in the tt invariant mass spectrum, the details (peak position and width) 
of course, depending on the compactification scale and the effective coupling K./M p i. 
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5.2 Single top quark production 



The interest in single-top production at hadron colliders comes from a number of reasons. It 
allows for studies of charged-current weak interactions of the top quark and for a direct extraction 
of the CKM-matrix element V t t,. Moreover, depending on the model under consideration, the cross 
section for single-top production acquires large corrections in BSM scenarios. At Tevatron first 
evidence for single-top production has been found only rather recently [140, 141]. 



Figure 16: Sample of Feynman diagrams for single top-production. Shown are the production in the s- 
channel (left), in the ^-initiated ^-channel (second from left), in the ^-initiated ^-channel (second from right) 
and in the &g-channel (right). 

The distinct channels for single-top production are displayed in Fig. [161 where samples of lead- 
ing order Feynman diagrams are shown. The ^-channel mode (Fig. |T6l left) proceeds through the 
production of an off-shell W-boson and subsequent decay (W — > tb). The initial state is propor- 
tional to the light flavor PDFs and the rate at LHC is relatively small. The f -channel exchange 
of a virtual W-boson in the boson-gluon fusion mode (Fig. [TH second from right) on the other 
hand is the dominant production mechanism, both at Tevatron and LHC. In the latter case, it is 
much enhanced due to the gluon PDF and subsequent splitting to a bb-pair (g — > bb), while the 
^-initiated /-channel process ub — > dt itself (Fig. [161 second from left) is suppressed by the nu- 
merically small bottom PDF. Of course, precise predictions at higher orders require a consistent 
matching of both processes, i.e. g — > ^-splitting in the hard parton scattering and in the evolu- 
tion [142]. For this main single-top production mode the NLO QCD corrections have also been 
subject of the latest addition to MC@NLO [143], so that NLO parton level calculations and show- 
ering are consistently combined. Finally, there is WY-production in the &g-channel (Fig. [161 right) 
being the second largest mode at LHC, but with negligible rates at Tevatron due to limited phase 
space. The £>g-channel is subleading in QCD (counting powers of the coupling constants) but again 
enhanced by the gluon luminosity. 

Sensitivity to BSM models appears in different manifestations for the various channels of 
single-top production. If one allows for anomalous couplings or flavor changing neutral currents, 
the f -channel contribution is altered significantly. On the other hand, the existence of charged "top- 
pions", excited Kaluza- Klein modes of the W-boson or, similarly, a W'-boson would have impact 
on the ^-channel. 
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5.3 Top quark plus jet production 



The final example of this Section is concerned with top quark plus jet production where, due to the 
high center-of-mass energy at LHC, also large statistics is expected. At the same time, the process 
tt+ jets is an important background to Higgs or supersymmetry searches, so that experimental 
search strategies need to impose kinematical cuts on the final state. At tree level, the cross sections 
have a really large scale dependence, which is why NLO QCD corrections are mandatory. For 
instance, the process tt + 2 jets entered Tab.[2]as a specific background to ttH, where the Higgs de- 
cays into a bb-pak. The NLO QCD corrections to the former would help to control the background 
uncertainty due to a heavy flavor mistag in a tt + 2 jets event. 
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Figure 17: The scale dependence of the LO and NLO cross sections for tt + 1 jet production at the 
Tevatron (left) and at LHC (right) with renormalization and factorization scales identified, fi r = /jf = /j 
(from Ref. [144]). 



On the way to this challenge (six-leg processes being currently at the edge of technology), the 
process tt + 1 jet was computed recently to NLO QCD in an impressive state-of-the-art calcula- 
tion [144]. Fig. [FT] displays the much improved scale dependence and shows that the perturbative 
corrections are moderate for the nominal scale choice ^i^m t . Clearly, it will be very interesting to 
see the NLO differential distributions for this reaction in the future and to compare them with LO 
predictions (e.g. from MadGraph) combined with parton showers. In this way, one can assess how 
well the NLO predictions for jet-observables are modeled by the underlying partonic processes 
such as gg — > ttg. 
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6 Higgs production at LHC 



Let us conclude this review with a brief discussion of the flagship measurement to be conducted at 
LHC. To start with, we illustrate in Fig.[T8]for the Standard Model Higgs the dominant production 
modes (left) and the branching ratios for the decay (right) as a function of the Higgs mass. The 
plotted values for the mass range up to M# = 1 TeV, which is generally considered an upper 
bound for the Standard Model Higgs due to triviality. A lower bound on the Higgs mass has been 
established from direct searches at LEP [145, 146], currently M# = 114.4 GeV. With the high 
luminosity and statistics of run-II the Tevatron experiments currently conduct an active search for 
the Higgs as well (see [147]). 

Focusing on hard QCD aspects, we will limit ourselves to the production part (i.e. Fig. [T8] 
on the left). Depending on the Higgs mass, the signatures from the various decay modes (gauge 
bosons, lepton pairs, quark pairs, etc.) define the experimental search strategy and, at the same 
time of course, the need to improve predictions for the competing Standard Model processes as 
discussed in Sec. |4l Unfortunately, this will not always be successful, as for example the fully 
hadronic modes will not be accessible for Higgs detection because of the huge QCD mult-jet 
background. However, due to the numerous channels we can hardly touch these aspects. A very 
extensive discussion of Higgs physics at colliders can be found e.g. in the recent Refs. [148, 149]. 




M H [GeV] M H [GeV] 

Figure 18: Left: the total cross section for Higgs production at LHC at NLO in QCD. Right: the branching 
ratios of Higgs boson decay (from Ref. [148]). 

As can seen from Fig. [T8]on the left, we have a clear hierarchy of channels and in Fig. [19] we 
display samples of Feynman diagrams for these production modes. Gluon fusion is induced via 
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Figure 19: Sample of Feynman diagrams for the various modes of Higgs production. Gluon fusion (left), 
weak vector-boson fusion (second from left), Higgs-strahlung (second from right) and ??H-channel (right). 



a heavy quark loop (Fig. [19] left). It has the largest rate for all values of the Higgs mass Mh due 
to the large top-Yukawa coupling of the Higgs boson and the large gluon luminosity. In limit of a 
heavy top (m t — > °°) one can describe the interaction by an effective ggH-vertex upon integrating 
out the heavy quark in the loop, which is a very good approximation also for finite m t . QCD 
corrections for numerous observables in this channel have been determined and we will high-light 
a few aspects below. 

Weak vector-boson fusion proceeds via qq — > qqH (Fig. [19] second from left) and is mediated 
by ^-channel gauge boson exchange. It has the second largest rate, being dominated mostly by 
the u, J-quark PDFs and is proportional to the WWH coupling. However, the signal identification 
for the three-body final state needs dedicated cuts on the final state jets (see e.g. [150, 151]). The 
characteristics of the latter are extremely important to discriminate the VBF signal from QCD 
backgrounds, for instance through forward jet-tagging and central jet-vetoing (see e.g. [152]). 
Higgs-strahlung in the channel qq — > W(Z)H (Fig. [19] second from right) makes up for the third 
largest rate and has the same couplings as vector boson fusion. The production mechanism requires 
a gauge boson from ^-annihilation so that the radiative corrections in QCD for Higgs-strahlung 
follow largely from the corresponding ones in the Drell-Yan process and have been determined to 
NNLO [153]. 

As another mode the associated Higgs production with heavy quarks has been discussed, for 
example pp — >■ ttH (Fig. [19] right), which is now known to NLO in QCD [154-157]. At LHC 
the process is driven by the gluon luminosity but the rate drops quickly for larger Higgs masses 
and the phase space becomes too small already for values of Mh — 180 GeV. In addition, for 
lower Higgs masses the final state from pp — > ttH has a large Standard Model background (see 
Tab. [2]) which will be difficult to suppress. The process pp — > bbH even has a slightly larger 
rate at LHC for Mh < 300 GeV, but the final state with ^-jets is overwhelmed by background. 
Generally, Higgs couplings to bottom quarks are more important in the extended Higgs sector of 
the MSSM (see e.g. [149]). Finally, the production of Higgs boson pairs {pp — > HH +X) at LHC 
still has rates accessible to measurements for low Higgs masses, however we will not discuss these 
processes further. Cross sections for the production of three or more Higgs bosons at LHC are too 
small [148]. 

Let us, for the rest of this Section focus on the role of higher order QCD corrections for the 
Higgs signal in gluon fusion. The total cross section for gluon fusion in the heavy top limit is 
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known exactly to the NNLO in QCD [158-160]. The radiative corrections yield a sizable ^-factor 
of two or more as can be seen in Fig .1201 (left) where we display the total cross section as a function 
of Mh in the low mass range. Of course, the origin of this large effect is well understood. It is due 
to soft gluon emission, which makes up numerically for the bulk of the perturbative corrections (see 
e.g. [161] and references therein). This fact has motivated the derivation of the complete soft N 3 LO 
corrections in [161] which are also plotted in Fig. [20] and illustrate nicely the property of apparent 
convergence of the perturbative expansion. Another indicator in this respect is, as often stressed so 
far, the stability under scale variation in a typical range, say Mh < /J < 2M#. This is shown in in 
Fig. [20] (right) where we display the total cross section as a function of the renormalization scale 
for a Higgs mass Mh = 120 GeV. As the scale dependent terms are completely predicted by the 
lower order terms, the curve denoted by N 3 LOapprox should be a very good approximation of the 
exact three-loop result with a residual uncertainty estimated of O ( 1 — 2%) only. 



o(pp -> H+X) [pb] 
Hr = m h 



N-LO 




N LO„ 




- 40 - 
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Figure 20: Perturbative expansion of the total cross section for Higgs production at the LHC. Shown are 
the dependence on the mass Mh and the renormalization scale /j r (from Ref. [161]). 



Comparison with experimental data of course requires selection cuts on differential distribu- 
tions. The latter are available for the gluon fusion channel including QCD predictions up to NNLO 
and allow to study bin-integrated distributions (e.g. for the Higgs rapidity) with subsequent Higgs 
decay in a variety of modes. Presently the NNLO corrections in gluon fusion have been combined 
with the decay modes H -> yy as well as H -> WW -> Mv and H -> ZZ -> 41 [28, 29, 162, 163]. 
It is very interesting to study these higher order corrections under the impact of kinematical cuts 
on the observed final state leptons, photons or the jet activity. Contrary to the findings for the total 
cross section, where higher order corrections amount to 100% or more, the effect of radiative cor- 
rections in distributions is strongly reduced by the selection cuts. We illustrate this for the Higgs 
rapidity distribution as calculated with the parton level Monte Carlo program HNNLO [29]. In 
Fig. EH on the left no cuts on the p t of additional jets are applied and the increase of the NNLO 
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over the NLO corrections amounts to approximately 20% with some dependence on the rapidity. 
In contrast, in Fig. |2H on the right all jets with p t > 40 GeV have been vetoed, a situation typi- 
cal when searching for the Higgs in the decay mode H — > WW to suppress the WW background 
from tt- production. As an upshot, the size of NNLO QCD radiative corrections is reduced to 5%. 
Moreover, as expected, the NNLO QCD corrections improve significantly the stability under scale 
variation. First steps towards an assessment of the WW-background from QCD at NNLO have 
recently been made [164]. Also we remark that for the Higgs rapidity distribution even the soft 
N 3 LO corrections have been obtained and simple analytical formulae is available [165]. 

In summary we conclude that the rates for the Standard Model Higgs at LHC are reliably pre- 
dicted by QCD. We have illustrated this for the gluon fusion channel where we have observed how 
higher order QCD predictions decrease the sensitivity to scale variations. There exists a residual 
uncertainty of the cross section of a few per cent due to the parton luminosity (see Sec. 13.41) . How- 
ever, the gluon PDF is well constrained in the kinematical range and HERA data can be evolved 
to NNLO accuracy. The study of differential distributions in particular with realistic experimental 
cuts is an active field of ongoing research for all main production modes (see Fig. [T9l) an we refer 
to the literature (see e.g. [7, 8, 62, 148]). 



7 Summary 

We have briefly reviewed the theoretical framework of QCD at hadron colliders. Precision predic- 
tions for hard scattering cross sections rely on a detailed knowledge of the parton content of the 
proton and of the rates for the corresponding partonic subprocess. We have given an overview of 
our current information on parton distributions including evolution to LHC energies. For various 
Standard Model cross sections such as W ± /Z-boson, Higgs boson or top quark production we have 
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reported on the present status of perturbative QCD predictions. Precision cross sections require 
the calculation of higher order corrections for both, the signal and the background with massive 
particles and jets. A lot of ongoing activity in this direction is concentrated on processes with 
multi-particle production and we have tried to give a snapshot of the technology involved, e.g. 
the idea on- shell recursions for scattering amplitudes. In summary, we have tried to convey the 
message that QCD theory is ready to meet the challenges of LHC. 

In this review, we have mostly omitted details of specific hadronic final states, e.g. jet algo- 
rithms, /?-quark (£»-jet) production or aspects of Z?-quark fragmentation as well as parton showers in 
Monte Carlo simulations. We have also left out any discussion of resummation approaches meant 
to improve fixed order perturbation theory, be it threshold logarithms of Sudakov type or \n(p t )- 
terms in transverse momentum. Finally, nothing has been said about the region of small-* and for- 
ward physics at large rapidities (e.g. for diffractive production of Higgs bosons). For all these re- 
maining aspects as well as a broader coverage, the interested reader is referred to [2,3,7,8,37,38,62] 
and the numerous references therein. 

QCD theory at hadron colliders is an extremely dynamical field at the moment, thus we expect 
that many specific issues will be improved or further clarified soon. The broad theoretical frame- 
work however, will certainly remain valid and hopefully, the present compilation has high-lighted 
the role of QCD in the era of LHC. 
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